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Torrefaction is a thermal treatment applied to improve biomass properties for long-term storage and 
thermochemical conversion of wood. The torrefaction kinetic behavior of two woody biomass samples 
(aspen and birch) and their individual constituents (cellulose, xylan -hemicellulose proxy- and lignin) 
was studied in detail. The materials were torrefied in a fixed-bed reactor up to 280 °C and were charac¬ 
terized after 15 min and 60 min treatments in terms of their energy density, energy yield, H, O and C con¬ 
tents, humidity uptake and FTIR vibrational spectra. Regardless of mildness (240 °C) or severity (280 °C) 
of torrefaction, hemicellulose was the most reactive. Compared to torrefied cellulose, lignin exhibited 
considerable activity and increased heating value under severe torrefaction (280 °C). Likewise, the 
increase in heating value of torrefied cellulose was less significant than that of lignin due to cellulose 
energy loss because of production of tar compounds. The release of volatiles versus time for cellulose, 
xylan and lignin torrefaction reactions measured using a thermogravimetric method was also modeled 
for each lump by means of the distributed activation energy model and a three-parameter log-normal 
distribution function. The model was found to fit the experimental data for xylan and lignin torrefaction 
while it showed slight under-prediction in the case of cellulose torrefaction. The distributed activation 
energy models developed for the individual lumps were combined into a meta-model to apprehend 
the torrefaction kinetics of actual (aspen and birch) woody biomass by assuming parallel and indepen¬ 
dent involvement of their three constituents present under various compositions. The model was found 
to represent with very good agreement the experimental data of aspen and birch torrefaction. 

© 2013 Elsevier B.V. All rights reserved. 


1. Introduction 

Environmental and economic constraints as well as resource 
scarcity are increasingly calling into question the classical fossil- 
energy driven development paradigm. In response, alternative 
solutions such as non-edible biomass are steadfastly climbing the 
ladder of the energy portfolio as potentially viable and perennial 
carbon neutral energy options [1,2], However, as a fuel lignocellu¬ 
losic biomass is disadvantaged in terms of O/C ratio, caloric value 
and moisture content, to name its main handicaps [3], In order 
to upgrade biomass to a higher quality fuel, thermal treatments 
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such as torrefaction are required. The goal behind torrefaction is 
to eliminate oxygen in the form of low caloric value torrefaction 
gases, e.g., CO 2 , water and some organic acids, by exposing biomass 
to temperatures as high as 200-300 °C [4-6], In spite of the fact 
that ca. 30 wt.% of biomass is lost during torrefaction the recovered 
solid may retain up to 90% of the initial biomass energy content [3], 
In comparison to the parent raw biomass, torrefied solid products 
should bear certain properties such as lower moisture content, 
lower O/C ratio, and higher energy density to make them suitable 
biofuels [4-6], 

Different torrefaction processes, e.g., dry torrefaction [4-6], wet 
torrefaction [7,8] and ionic-liquid assisted torrefaction [9] have 
been proposed to leverage the above biomass properties. Dry 
torrefaction consists of slow heating of biomass in an inert 
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atmosphere from 200 °C to 300 °C [4-6,10], During this treatment, 
the three main constituents of biomass, i.e„ hemicellulose, cellu¬ 
lose and lignin [11] decompose at different rates. Hemicellulose 
is the most sensitive lump to be affected by dry torrefaction unlike 
cellulose whose decomposition is the rate limiting step of torrefac¬ 
tion [12], The rate of cellulose decomposition is stimulated by ele¬ 
vated torrefaction temperatures, though temperatures in excess of 
300 °C are not recommended due to the release of heavy tars 
which tendency is to reduce biomass energy density [12], Impreg¬ 
nation of biomass with room-temperature ionic liquids enhanced 
the rate of cellulose torrefaction at temperatures lower than 
300 °C [9], However, despite advantages showcased by ionic liq¬ 
uids, the solid products resulting from ionic-liquid torrefaction 
treatments still require further improvements such as energy yield 
which was found to be lesser compared to the dry-torrefaction 
analogs [9], In addition, cost of ionic liquids should be factored in 
if such advanced approaches are to be sought for commercial 
applications of biomass torrefaction [13], Wet torrefaction or 
hydrothermal pretreatment of biomass performed in a pressurized 
liquid water up to 5 MPa and 260 °C produces solids with energy 
density greater than that obtained from dry torrefaction [7,8], Like¬ 
wise, elevated pressures tend to reflect in increased operational 
costs also. Therefore, it is likely that dry torrefaction, because of 
its rustic character and low cost, will still encompass the majority 
of lignocellulosic biomass torrefaction units [14], 

Several kinetic models have been proposed for biomass torre¬ 
faction in dry inert atmosphere which are mostly of pseudo-mech¬ 
anistic nature where one or several pseudo-reactions are assumed 
for modeling biomass torrefaction [12,15,16], We are of the opin¬ 
ion that the wide range of heterogeneous reactions taking place 
during torrefaction must call for more appropriate mathematical 
representations. These latter must be able to account for the dis¬ 
tributed nature of the reactions across the biomass structure, e.g., 
catalytic effects of the alkali metals, and the corresponding wide 
range of activation energies [17] which are not reflected in current 
models. We propose to evaluate the potential of the distributed 
activation energy model (DAEM) which is viewed as one of the 
most accurate approaches for modeling biomass pyrolysis [17], 
The model assumes an infinite number of irreversible nth-order 
parallel reactions with different activation energies occurring 
simultaneously and where the difference in activation energies is 
represented by a distribution function. The DAEM has been widely 
used to analyze complex reactions such as pyrolysis of various 
ranks of coal and other fuels [18-21], Yet, and to the best of the 
authors’ knowledge, the merit of using DAEM has not been 
assessed for representing biomass torrefaction kinetics. 

This work, therefore, aims at studying the dry torrefaction prop¬ 
erties and kinetics of the main biomass constituents, i.e., xylan 
(hemicellulose proxy), cellulose and lignin, as well as two different 
hardwoods, i.e., aspen and birch. Characterization experiments 
were designed such that they evaluate a number of properties 
which are of paramount importance in the field of biomass torre¬ 
faction and provide a rather consolidated study. To do so, the main 
biomass constituents, birch and aspen were torrefied under mild 
and severe torrefaction conditions and the contribution of each 
biomass constituents on the properties of torrefied biomass was 
analyzed by obtaining their heating value and energy yield. Major 
emphasis was put on lignin which is considered as an inert com¬ 
pound in some existing literature [22,23] and thus its decomposi¬ 
tion rate, heating value and energy density of torrefied lignin was 
compared with that of cellulose. The FTIR analysis performed on 
the raw and torrefied samples also helped gaining an understand¬ 
ing on the mechanism of cellulose, hemicellulose and lignin torre¬ 
faction. Moreover, the humidity uptake of each biomass 
constituents subjected to the mild and severe torrefaction condi¬ 
tions was studied to investigate their impact on equilibrium mois¬ 


ture content of torrefied birch and aspen. Finally, a DAEM kinetic 
model for the torrefaction of each constituent, separately, and for 
both hardwood samples is formulated. To achieve this, separate 
DAEM models are developed for each one of hemicellulose, cellu¬ 
lose and lignin lumps by monitoring their thermal decomposition 
rates versus time in a torrefaction environment using a thermo- 
gravimetric analyzer (TG). The quality of the models is judged both 
by their predictive power for actual wood samples using recom¬ 
bined kinetics from the constitutive wood lumps and actual torre¬ 
faction kinetics of aspen and birch samples. 

2. Experimental 

Isothermal torrefaction studies were carried out using trem¬ 
bling aspen and birch, two abundant woody biomass samples from 
Quebec forests (Canada), as well as powdered samples of cellulose 
(Avicel®PH-101), xylan (X0502) and lignin (37,059-lignin, alkali), 
all three purchased from Sigma Aldrich. The proximate and ulti¬ 
mate analyses of the studied samples are summarized, respec¬ 
tively, in Tables 1 and 2. The birch and aspen samples were first 
grinded and then sieved to isolate particles finer than 500 pm. 
The samples were then torrefied in a fixed-bed torrefaction reactor. 
The bench scale biomass torrefaction unit consisted of a stainless 
steel (SS) fixed-bed reactor (ID. 34 mm, length 70 cm) placed inside 
an induction coil. For each one of the experiments, 300 mg of bio¬ 
mass was held atop of a steel grid located in the middle of a stain¬ 
less steel tube to act as a support for the biomass bed. The reactor 
was heated through a 30 cm long coil connected to a 18 kW induc¬ 
tion power supply for heating the stainless steel tube via an AC 
magnetic field then heating in turn the biomass layer through radi¬ 
ation and conduction. The biomass bed temperature was raised 
from the ambient to the target temperature (240 °C and 280 °C) 
in less than 10 s and then the samples were kept at the final tem¬ 
perature for specified time durations, 15 and 60 min, respectively. 
These torrefaction temperatures, i.e., 240 °C and 280 °C were rep¬ 
resentatives of mild and sever torrefaction conditions, respectively 
[24], The SS tube temperature was controlled by a fiber optic 
pyrometer connected to a temperature controller to meter the 
induction-heater power. Temperature was also measured through 
a K type thermocouple inserted in the center of the biomass bed. 
The difference in temperature between the bed and SS tube surface 
never exceeded 5 °C in steady-state conditions. Volatiles forming 
in the reactor were directed towards a cold trap by sweeping the 
bed with a continuous N 2 flow (200 N mL/min). 

The solid products obtained after torrefaction were further ana¬ 
lyzed to quantify their properties after thermal treatment. The so¬ 
lid products were weighted to obtain their mass yield (Y M = %mass 
solid product/mass feed) after torrefaction. The C, H, N, S contents 
of the samples were also measured using an elemental analyzer 
(Fisons EA 1108CHNS). All measurements were repeated twice 
and their averages were reported. The higher heating value 
(HHV) of the samples before and after torrefaction were deter¬ 
mined using a relationship obtained by Channiwalaet al. [25] based 
on carbon, hydrogen, oxygen, sulfur, nitrogen and ash content of 
biomass. The model developed in their work was tested for 225 


Table 1 

Proximate analysis in mass% of cellulose, xylan, lignin, birch and aspen wood. 


Humidity Volatile (dry Fixed carbon (dry Ash (dry 

basis) basis) basis) 


Cellulose 3.4 
Xylan 8.5 

Lignin 3.1 

Birch 5.5 


93.1 

74.2 

49.3 

89.4 

87.5 


5.9 


47.5 

9.8 

12 
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te analysis, O/C and H/C m 


ild, high heating value, energy yield and energy density of as re 


Mass yield (%) Energy yield HHV (kj/g) Energy density ra 


42.8 6.4 0 


0 6.0 0.1 


50.8 0.89 


53.4 0.99 


38.9 0.44 


6.0 0.90 2.29 27.2 0.32 


5.5 0.60 2.12 22.2 


46.7 0.74 


43.9 0.62 


43.0 0.59 


samples and was shown to be applicable for estimation of the high 15 and 180 min while samples’ mass was adjusted so that weight 
heating values of torrefied biomass [26,27]: losses during torrefaction never dived below 30 wt.%. 


HHV (kj/g) = 0.3491 x C + 1.1783 x H + 0.1005 x S - 0.1034 
x 0-0.0151 xN-0.0211 x Ash 


With C, H, S, 0, N and Ash representing carbon, hydrogen, sulfur, 
oxygen, nitrogen and ash mass percentages, respectively. Using the 
solid product HHV, the energy density (e = HHV solid product/HHV 
feed) and energy yield (Y e = e f M ) of the solids were also obtained. 

To investigate the effect of torrefaction on the samples hydro- 
phobicity, equilibrium moisture content (EMC) analysis was also 
performed on the torrefied samples. EMC of raw and torrefied bio¬ 
mass samples were obtained by putting 300 mg samples in a va¬ 
por-saturated environment controlled at 25 °C. All the samples 
were oven-dried at 110 °C overnight to remove the pre-existing 
humidity in the samples before EMC analysis. The samples weight 
gain was measured after ten days which was estimated to be 
enough for the samples to reach moisture equilibrium uptake. Fou¬ 
rier transform infrared spectroscopy (FT1R) (Nicolet Magna 850 
spectrometer, Thermo Scientific, Madison, WI) was also performed 
on the as-received and torrefied cellulose, xylan and lignin to study 
changes in their functional groups during torrefaction. The FTIR 
was equipped with high temperature Golden Gate ATR accessory. 

For the derivation of the kinetic model for biomass torrefaction, 
isothermal torrefaction kinetics of all five substrates were per¬ 
formed in a thermogravimetric analyzer (TG) (Perkin Elmer Lab 
System Diamond TG-DTA) used to monitor the weight loss of the 
samples during torrefaction. TG is appropriate for developing ki¬ 
netic models for steady-state (isothermal) or low-heating rate bio¬ 
mass torrefaction [28], For this purpose, the samples were placed 
in an alumina ceramic holder forming a small bed with initial mass 
of 5-8 mg of woody biomass, lignin, cellulose, or xylan. The sam¬ 
ples were first dried in TG at 120 °C for 15 min and then heated 
up to the target temperature (240-290 °C) with the maximum 
achievable heating rate (-200 °C/min) under N 2 flow (200 N mL/ 
min). The isothermal torrefaction of the samples was performed 
by maintaining the samples at the target temperature between 


3. Kinetic modeling 


The total amount of volatiles released during torrefaction of 
xylan, lignin and cellulose were modeled separately using DEAM 
from the mass loss registered via TG. In this model, the decompo¬ 
sition of xylan, lignin and cellulose during torrefaction is consid¬ 
ered to result from an infinite number of n th-order parallel 
reactions. For a particular ith reaction decomposition rate, one has: 


m/vi) 
a - w 


= kdt 


(1) 


where “i” may contribute to the decomposition of xylan, lignin and 
cellulose. In Eq. (1), Vi refers to the volatile contribution from reac¬ 
tion “i” during torrefaction of each biomass model compound and 
V* denotes the maximum volatile contribution to be released due 
to reaction “i”. The rate constant of reaction “i” is expressed follow¬ 
ing an Arrhenius form as: 

ki = A,exp( -Ej/RT) (2) 

In which R is the universal gas constant, A t is the frequency fac¬ 
tor, E; is the activation energy and T is the torrefaction tempera¬ 
ture. Considering that the number of parallel reactions during the 
torrefaction of each compound is infinite, the maximum volatile 
contribution released by each reaction (V?) can be considered as 
a fraction of the total volatiles (V*) with activation energies be¬ 
tween E, and E,- + <5E ; , i.e., V* = V*/(E,■)<■)£;). Consequently, the total 
volatiles evolved up to time instant t during torrefaction of each 
compounds can be expressed as [29]: 


^ = 1 - Y = /°V( £ ) [!-(!- n) J‘ Aexp(-E/RT)dt\ ' d£ (3) 

In Eq. (3), J[E) is the distribution of activation energies 
(JiT /(£)d£ = 1), V is the amount of volatiles evolved up to time t, 
V* is the maximum volatile content of each compound and V is 
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the residual compound volatile which is not released. In the special 
case where the order of the reaction is equal to one, Eq. (3) reduces 
to: 

£ = i - r = [me- W-WdE (4) 

Several probability density functions such as Gaussian [18,29], 
Weibull [30,31] and Gamma [32] distributions have been tested 
for biomass pyrolysis. In the present study, the three -parameter 
log-normal distribution was used to model the distribution of acti¬ 
vation energies of cellulose, xylan and lignin during torrefaction. 

f ° r £ > E " (5) 

f(E) = 0 for 0 <E<Eq 

Substituting/[£) in Eqs. (3) and (4) by a log-normal distribution 
function leads to a double integral for which there is no analytical 
solution. Alternatively, numerical methods such as recursive adap¬ 
tive Simpson quadrature technique available in MATLAB were 
used. There are five unknown parameters in the torrefaction DAEM 
for cellulose, lignin and xylan, i.e., /r, i5 and E 0 being the unknowns 
in the log-normal distribution and n and A associated with the or¬ 
der and frequency factor. These torrefaction kinetic parameters 
were obtained using a simulated annealing optimization tech¬ 
nique. This technique is a derivative-free optimization method rec¬ 
ommended for finding global optima for DAEM [33], 

The DAEM developed for the isothermal torrefaction of cellu¬ 
lose: xylan and lignin could be employed to develop a kinetic mod¬ 
el for the torrefaction of aspen and birch. It is assumed that the 
aspen and birch torrefaction is the result of torrefaction of their 
constitutive cellulose, xylan and lignin in a way that there is no 
interaction between components. Consequently, the following 
equation can be written to relate the torrefaction of aspen and 
birch at time t to the torrefaction of their main constituents. 

(vvv-)^~|^(v/n (6) 

In Eq. (6), q represents the contribution of each compound to 
the total volatiles release. Using the simulated annealing optimiza¬ 
tion technique, the optimum values of q which could fit the three 
component model to the experimental torrefaction data of aspen 
and birch could be obtained. 


4. Results and discussion 

4.1. Biomass torrefaction 

Lignin torrefaction. Lignin is a stable biopolymer which confers 
rigidity and strength to cellulose in the plant cell wall [34], It con¬ 
sists of different proportions of guaiacyl (G), p-hydroxyphenyl (H) 
and syringyl (S) units which are inter-bonded irregularly through 
different types of linkages; the so-called p0 4 linkage being the 
most frequent one [35-37], Lignin resistance to decomposition 
depends on the type of linkages between its building units [36], 
The aryl ether linkages (C—O—C) have low disassociation energy 
while the biphenyl type linkages are more resistant to breakage 
[36], Lignin decomposition takes place at temperatures as low as 
150 °C, earlier than xylan and cellulose though at quite a tepid rate 
(Fig. 1). Lignin decomposition rate increases slowly with tempera¬ 
ture until 300 °C and reflects in a slight decrease in solids mass 
yield by increasing the torrefaction temperature (Table 2). For in¬ 
stance, lignin torrefaction at 240 °C and 280 °C led, respectively, 
to 89.9% and 87.3% of dry mass yield. The lignin torrefaction is 
slightly exothermic as can be observed in Fig. 1 . 



cellulose, xylan, lignin, birch and aspen wood decomposed in thermogravimetric 
analyzer at 5 °C/min heating rate until 300 °C. 


Torrefaction reduces the O/C and H/C atomic ratios in lignin (Ta¬ 
ble 2). As an example, these ratios are reduced from H/C =1.12 and 
O/C = 0.32 in raw lignin to H/C = 0.93 and O/C = 0.23 for lignin torr¬ 
efied at 280 °C. Also, these ratios are sensitive to torrefaction tem¬ 
perature such that the lignin carbon content increases with 
increasing the torrefaction temperature, whereas those of oxygen 
and hydrogen decrease. Reduction in O/C and H/C ratios can also 
be expected from the analysis of lignin decomposition volatiles 
at temperatures lower than 300 °C. The release of H 2 0, C0 2 , form¬ 
aldehyde and formic acid with higher O/C and H/C ratios than raw 
lignin, as a result of the cleavage of OH and carbon groups in the 
alkyl side chains, explains indeed these ratios’ trends in torrefied 
lignin [36]. The energy density of torrefied lignin at 240 °C and 
280 °C increases by 7% and 9% (on a dry basis), respectively, due 
to the reduction in the oxygen content of lignin (Table 2). Consid¬ 
ering that almost 3% of raw lignin consists of water, the energy 
density of torrefied lignin improves by ca. 10% with respect to that 
of untreated lignin. Lignin mass loss amounts to 13-17% due to 



Wavenumber (cm- 1 ) 


Fig. 2. FTIR spectra of as-received and torrefied cellulose, xylan and lignin samples. 
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torrefaction treatment while almost 95% of its energy is preserved 
in the torrefied lignin product (Table 2). According to the van Krev- 
elen diagram, raw lignin shares the fuel characteristics of low rank 
coals such as lignite or sub-bituminous coal [3], 

Lignin FTIR spectra before and after torrefaction at 240 °C and 
280 °C are shown in Fig. 2. The aromatic nature of lignin due to 
the guaiacyl (G), p-hydroxyphenyl (H) and syringyl (S) units man¬ 
ifests in the vibrational range of 1400-1600 cm -1 due to the C—C 
bonds in the aromatic structure [35]. Likewise, both aliphatic and 
phenolic hydroxyl groups, which prompt hydrogen bonding in 
the lignin polymeric structure, exhibit a broad hydroxyl peak in 
the 3200-3600 cm -1 range [38], Though torrefaction promotes lig¬ 
nin dehydroxylation and thus a reduction in the extent of hydro¬ 
gen bonding, such changes are barely visible in terms of 
broadness and relative of OH vibrational peaks in the FTIR spectra 
of the torrefied lignin (Fig. 2). Lignin torrefaction leads to a de¬ 
crease in the peak intensities near 1035, 1205 and 1265 cm -1 
which are assigned to C—O bond deformations in methoxyl groups, 
C—O bond stretching in phenolic hydroxyl groups and C—O bond 
stretching in the ether aromatics, respectively [35]. 

Torrefaction reduces the equilibrium moisture content of lignin 
as well. The major mechanism for water adsorption on lignin is 
water physisorption via hydrogen bonds between water and the 
free hydroxyl groups in the lignin subunits [39], Torrefaction also 
leads to a partial destruction of the free hydroxyl groups in the lig¬ 
nin. This translates in almost 40% reduction in equilibrium mois¬ 
ture uptake of treated lignin as shown in Fig. 3 Increasing the 
torrefaction temperature could extend the lignin dehydroxylation 
causing less water uptake by lignin. 

Cellulose torrefaction. Cellulose consists of linear chains of sev¬ 
eral hundreds to thousands of jS(l -»4) linked D-glucose units 
[40,41], Cellulose decomposition commences around 260 °C with 
a rate that increases with temperature (Fig. 1 ). Similar to lignin, 
the torrefaction of cellulose is slightly exothermic. Torrefaction of 
cellulose at mild and severe conditions led to 1.9% and 9.9% mass 
reduction as can be seen in Table 2. The lower mass yield of torr¬ 
efied cellulose at mild and severe conditions in comparison with 
lignin shows its lower rate of decomposition than lignin in the tor- 
refaction temperature range. The mass reduction of mild condition 
torrefied cellulose is considerably lower than that of torrefied 
lignin at the same torrefaction conditions (1.9% and 10.1%, respec¬ 
tively). Although cellulose decomposition at mild condition torre¬ 


faction was low, a reduction in H/C and O/C ratios in cellulose 
was observed (Table 2). The mild torrefaction led to the reduction 
in the H/C and O/C ratios of cellulose from H/C =1.81 and 
O/C = 0.89 in raw cellulose to H/C = 1.72 and O/C = 0.86 in the mild 
condition torrefied cellulose. The decrease in the H/C and O/C 
ratios of mild condition torrefied cellulose could be due to the in- 
tra-molecular elimination of water from the torrefied cellulose. 
Due to the reduction in the O/C and H/C of torrefied cellulose, 
the energy density of mild condition torrefied cellulose was in¬ 
creased by 1%. In the case of severe condition torrefied cellulose, 
the mass yield of torrefied cellulose was more comparable with 
that of torrefied lignin in the same conditions (90.1% and 87.3%, 
respectively, Table 2). However, the increase in carbon content of 
torrefied lignin at the same conditions is quite higher than that 
of cellulose (10.8% and 3.7%, respectively). The same trend was also 
observed for the oxygen and hydrogen content of torrefied cellu¬ 
lose and lignin at severe torrefaction condition. The reduction in 
oxygen and hydrogen content of severe condition torrefied cellu¬ 
lose were 2.75% and 3.13% while oxygen and hydrogen content 
of lignin reduced by 21.3% and 8.3%, respectively. These results 
indicate that the severe condition torrefaction of cellulose leads 
to the loss of some carbon in the form of tar. The release of tar 
could be the main reason for the small increase in the energy den¬ 
sity of severe condition torrefied cellulose than lignin. Although 
some previous torrefaction studies [22,23] indicated that lignin is 
relatively inert under torrefaction conditions, these results reveal 
some activity which might not be neglected in comparison to that 
of cellulose at equal temperature and treatment time. 

FTIR spectra of the as-received cellulose and the torrefied cellu¬ 
lose at 240 °C and 280 °C are shown in Fig. 2. The changes in the 
FTIR spectrum of torrefied cellulose at 280 °C (severe treatment) 
is more pronounced than that for 240 °C (mild treatment). How¬ 
ever, the cellulose main vibrational features are still perceptible 
after torrefaction. Cellulose torrefaction at 280 °C leads to partial 
dehydroxylation reflecting in decreased intensities of the vibra¬ 
tional peak at 3200-3400 cm 1 (O—H stretching) and 2850- 
2950 cm -1 (C—H stretching). The FTIR spectrum of torrefied cellu¬ 
lose at 280 °C pinpoints to the inception of two new peaks at 
1620 cm 1 and 1705 cm 1 attributed, respectively, to unsaturated 
carbon-carbon (C=C) and carbonyl (C=0) stretching bonds. The 
dominant mechanism for cellulose decomposition at torrefaction 
temperatures below 300 °C occurs via cellulose dehydration [42], 



Fig. 3. Humidity uptake of the 


id torrefied cellulose, xylan, lignin, birch 


ig 5 intervals of 48 h. 
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Intramolecular elimination of water from pyranose-ring C 2 and C 3 
carbons leads to unsaturated pyranose rings (enol-anhydrocellu- 
lose) recognizable from the incipient vibrationals peaks at 
1705 cm -1 as well as formation of ketone species (keto-anhydro- 
cellulose) as revealed by 1620 cm 1 peak [42], However, at torre- 
faction temperatures over 250 °C, in addition to dehydration 
reactions, cellulose degradation could also cause the formation of 
some levoglucosan (tar), CO, C0 2 and some acid hydrocarbons 
[43], Cellulose torrefaction favors loss of polar hydroxyl groups 
leading to thus to lesser water adsorption via hydrogen bonding. 
This is confirmed from Fig. 3 where nearly 30% reduction in equi¬ 
librium moisture content of torrefied cellulose was noted in com¬ 
parison to moisture uptake of raw cellulose. 

Hemicellulose torrefaction. After cellulose, hemicellulose is the 
most abundant polysaccharide group in plant [44], Depending on 
its source, it could contain different sugar monomers such as glu¬ 
cose, xylose, mannose, galactose, rhamnose and arabinose 
[45,46], According to the type of sugar residues, hemicellulose spe¬ 
cies are grouped into four classes such as xylans, xyloglucans, and 
mannose and mixed-linkage /i-glucans [47], Each class could be 
also divided into several subgroups depending on the type and 
chemical nature of the side chains in their structure. Hardwood 
hemicellulose is mainly composed of xylan which has a linear 
polymeric skeleton form constituted of 1 -4 liked /j-xylano pyran- 
osyl units which could be substituted by O-methyl-D-glucuronic 
acid residue [48], The acetyl groups could also substitute the side 
chains in the xylanopyranosyl units [48,49], 

Xylan decomposition initiates at temperatures around 200 °C. 
This temperature is lower than that of cellulose but higher than 
that of lignin (Fig. 1). The lack of crystallinity and lower polymer¬ 
ization degree of hemicellulose as compared to cellulose could lead 
to less thermal and chemical stability for hemicellulose [49], The 
rate of xylan decomposition increases with temperature up to a 
maximum near 280 °C (Fig. 1). Therefore, temperatures lower than 
300 °C are suitable for xylan torrefaction. Xylan torrefaction under 
mild and severe conditions also led to mass losses, respectively, 
45.5% and 58.8% (Table 2). These levels are considerably higher 
than those of lignin and cellulose for the same torrefaction condi¬ 
tions. It can be concluded that hemicellulose is the component 
most sensitive to torrefaction treatments. Xylan energy losses 
due to mild and severe torrefaction represent some 28.9% and 
44.1%, respectively (Table 2). These loss percentages are lower than 
the corresponding mass loss percentages in the same conditions. 
The high oxygen concentration in xylan lowers its energy density 
and thus its higher heating value. However, torrefaction reduces 
oxygen concentration in xylan and eventually increases its energy 
density. The loss of oxygen from xylan upon torrefaction could also 
be the reason for its slight exothermic decomposition (Fig. 1). The 
H/C and O/C ratios decreased, respectively, from 1.82 and 0.99 in 
the as-received xylan to 0.92 and 0.51 (mild treatment) and to 
0.81 and 0.44 (severe treatment). This was mirrored by an increase 
of 1.3% and 1.4% of energy density of torrefied xylan, respectively, 
for the mild and severe treatments. 

The FTIR spectra of as-received and torrefied xylan are shown in 
Fig. 2. It can be seen that xylan torrefaction at mild conditions 
could not complete the transformation xylan chemical structure. 
Xylan decomposition at temperatures lower than 250 °C is known 
to be preceded by the release of water, CO, C0 2 , methanol, alde¬ 
hyde type and acetic products [48], The peak at 1690 cm 1 as¬ 
signed to carbonyl (C=0) functional groups in the torrefied xylan 
is a marker to follow the progress of xylan decomposition at 
240 °C. However, xylan torrefaction at 280 °C (severe conditions) 
leads to major changes in xylan FTIR spectrum (Fig. 2). Decompo¬ 
sition of xylan at 280 °C is believed to lead to cleavage of glycosidic 
residues and release of furfurals, acids, water and other gases [48], 
Peak intensities in the ranges 950-1200 cm -1 (simple C—O bond 


stretching) and 3200-3600 cm _1 (simple O—H bond stretching) 
decrease considerably upon treatment indicating a loss of these 
groups. This is unlike the intensity of C=0 bond (1690 cm -1 ) which 
increased due to torrefaction at 280 °C. Likewise, xylan carboniza¬ 
tion during torrefaction and loss of hydrogen and oxygen led to the 
formation of aromatics (peaks at 1570 cm 1 ) as also confirmed by 
the spectral upward drift towards the large-wave number region 
[50], 

Fig. 3 shows that the moisture uptake by as-received xylan is 
considerably larger than that of cellulose and lignin. In addition 
to the presence of polar hydroxyl and carboxyl groups, the amor¬ 
phous structure of xylan enhances the adsorption of water in com¬ 
parison with cellulose. Upon torrefaction treatment, 
dehydroxylation and decarboxylation of xylan could lead to a con¬ 
siderable decrease in equilibrium moisture content. The equilib¬ 
rium moisture content of torrefied xylan at 240 °C and 280 °C is 
38% and 35% less than that of non-treated xylan. 

Birch and aspen torrefaction. Thermal decomposition of aspen 
and birch, as shown in Fig. 1, started at a temperature as low as 
180 °C likely via lignin decomposition. As discussed earlier, lignin 
decomposition took place at temperatures as low as 150 °C. The 
rate of aspen and birch decomposition slightly increased from ca. 
200 °C after decomposition of hemicellulose in birch and aspen 
was initiated, Fig. 1. Torrefaction of aspen and birch also exhibits 
a similar exothermic behavior as expected from the individual exo¬ 
thermic behaviors of their constituents. 

Results of elemental analysis as well as mass and energy yields 
of severe and mild torrefaction of aspen and birch are summarized 
in Table 2. The mild torrefaction of birch and aspen woods led to 
the respective mass yields of 86.2% and 86.4%. Increasing the torre¬ 
faction temperature also decreased the mass yield of torrefied 
birch and aspen to the respective values of 74.5% and 78.5%. The 
considerable decrease in mass yield of torrefied birch and aspen 
by increasing the torrefaction temperature was expected as the 
main constituents of birch and aspen such as cellulose, hemicellu¬ 
lose and lignin also showed the same trends due to torrefaction. As 
can be observed in this table, the H/C and O/C ratios of torrefied 
birch and aspen reduced substantially due to the torrefaction at 
mild and sever conditions. These reductions in H/C and O/C ratios 
of mild condition torrefied birch and aspen led to 7% and 6% in¬ 
crease in their energy density ratio, respectively. Although birch 
and aspen are only losing around 14% of their mass, the increase 
in their energy density is significant. At severe torrefaction condi¬ 
tions, the weight loss of birch and aspen increased to 25.5% and 
21.5%, respectively, but their energy density slightly increased 
about 9% and 7%, respectively, with respect to the non-treated 
wood samples. This trend could be explained due to the presence 
of high amounts of cellulose in the structure of birch and aspen 
woods. Almost 58.9% and 51.2% of birch and aspen wood are com¬ 
posed of cellulose [51,52], As observed during torrefaction of cellu¬ 
lose; the severe condition torrefaction of cellulose leads to release 
of some tar compounds such as levoglucosan [43] due to its rapid 
decomposition causing energy loss during torrefaction. Conse¬ 
quently, although the mass yields of severe condition torrefied 
birch and aspen woods are much lower than those of mild condi¬ 
tion torrefied birch and aspen, the difference in their energy den¬ 
sity is not that noticeable. 

The equilibrium moisture content of torrefied aspen and birch is 
reduced by 30% due to torrefaction, Fig. 3. This trend was expected 
from investigation of the individual biomass constituents. The loss 
of some biomass polar groups such as hydroxyl and carboxyl 
groups during torrefaction is the main reason for the observed 
behavior. The equilibrium moisture content of the as-received as¬ 
pen and birch was 30% and 26%, respectively. Torrefaction reduced 
its level to 17% for both wood samples. Torrefaction as a strategy to 
improve resistance to fungal colonization correlates with the 
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reduction of equilibrium moisture content of the treated wood. 
Since fungal growth on lignocellulosic materials is strongly corre¬ 
lated to water uptake [53], A significant wood decomposition 
was reported to occur above the wood fiber saturation point at 
25-35% moisture content [54], It should be mentioned that limited 
biomass decomposition due to fungal growth may also occur at 
moisture levels as low as 20% [55]. Certain types of fungi such as 
Aspergillus and Penicillium could grow at biomass moisture con¬ 
tents below 15% [56], Hence, attenuation of microbiological activ¬ 
ity in torrefied wood could facilitate its long-time storage. 

4.2. Kinetic modeling of torrefaction 

In the following section, first a DAEM was developed for the tor- 
refaction of the individual components, cellulose, xylan and lignin. 
The DAEM parameters for each compound was optimized so that 
the model fits the TG mass loss profile for each compound versus 
time at three different torrefaction temperatures, 240 °C, 260 °C 
and 280 °C. The DAEM optimal parameters for xylan, lignin and 
cellulose as gathered in Table 3. Furthermore, the DAEM was tested 
for predicting the TG mass loss profiles of torrefied xylan, lignin 
and cellulose versus time measured at 250 °C, 270 °C and 290 °C. 
As shown in Figs. 4-6, the DEAM fits well the instantaneous of Y 
/V* obtained by TG experiments in the temperature range of 
240-290 °C, respectively, for xylan (r 2 = 0.99), lignin (r 2 = 0.98) 
and cellulose (r 2 = 0.94). In the case of cellulose (Fig. 6), the DAEM 
model slightly under-predicts the instantaneous Y\Y profiles at 
250 °C and 270 °C while it shows small over-prediction at 280 °C 
and 290 °C. The optimum reaction orders for xylan, lignin and cel¬ 
lulose torrefaction reactions were 3, 2 and 1, respectively. 

Fig. 7 shows the activation energy distribution of cellulose, xy¬ 
lan and lignin torrefaction reactions obtained using the log-normal 
distribution function. As can be seen, activation energy distribu¬ 
tions of cellulose, xylan and lignin torrefaction reactions exhibit 
their respective maxima at 207 kj/mol, 132 kj/mol and 196 kj/ 
mol. In the case of cellulose and xylan, the activation energy distri¬ 
butions of the isothermal torrefaction reactions lie within the nar¬ 
row ranges of 200-220 kj/mol and 125-145 kj/mol, respectively, 
showing that the parallel reactions involved in the cellulose and 
xylan torrefaction have close activation energies. Narrow ranges 
of activation energy distribution for non-isothermal pyrolysis of 
cellulose and hemicellulose have been already reported. Sonobe 
and Worasuwannarak [17] derived the distribution of activation 
energy for non-isothermal cellulose pyrolysis at temperatures as 
high as 600 °C and a narrow range of activation energy (180— 
190 kj/mol) was reported with a maximum at 185 kj/mol. Gaspar- 
ovic et al. [57] also obtained an activation energy distribution for 
non-isothermal hemicellulose pyrolysis at temperatures as high 
as 500 °C using Gaussian distribution functions. They stated that 
the activation energy of hemicellulose pyrolysis is distributed in 
a range of 124-142 kj/mol with a maximum at 132.9 kj/mol. The 
narrow range of activation energy distribution for cellulose and xy¬ 
lan parallel torrefaction reactions suggests that a distributed acti¬ 
vation energy model for xylan and cellulose torrefaction might 
be regarded as equivalent to that in which only one single reaction 
represents the overall cellulose and xylan torrefaction. 



Time (min) 


Fig. 4. TG derived profile of volatile release versus time for xylan torrefaction and 
DAEM simulated profiles of xylan torrefaction. 



Time (min) 


Fig. 5. TG derived profile of volatile release versus time for lignin torrefaction and 
the DAEM simulated profiles of lignin torrefaction. 

In contrast to cellulose and xylan, the activation energy for the 
lignin torrefaction reactions is distributed widely in the range of 
116-384 kj/mol. This is in accordance with the wide activation en¬ 
ergy distribution of Alcell lignin pyrolysis (129-361 kj/mol) re¬ 
ported by Ferdous et al. [58] Mani et al. [33] also reported a wide 


Table 3 

Optimum values of parameters for the torrefaction kinetic models of cellulose, xylan, lignin, birch and aspen. 




E 0 (kj/mol) g (kj/mol) S (kj/mol) n 


89.8 3.76 

56.3 5.01 

190.8 2.82 


0.05 


0.15 


1.815 

1.516 

3.804 


10 16 

10 17 


0.534 0.235 0.231 

0.584 0.212 0.204 


0.99 

0.98 

0.94 
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0.97 
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Fig. 7. Distributions of activation energies for xylan, cellulose and lignin. 
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Fig. 9. TG derived profile of volatile release versus time for aspen torrefaction and 
the simulated trends using optimized c, values. 


range of activation energy for lignin pyrolysis (starting at around 
lOOkJ/mol) with maximum in the range of 158-170 kj/mol. The 
broad activation energy revealed for lignin torrefaction could be 
due to the presence of various functional oxygen groups in the 
structure of lignin with different thermal stabilities leading to par¬ 
allel reactions with broad activation energies. In addition, the pres¬ 
ence of alkali and alkaline metals found in the structure of lignin 
could also catalyse some decomposition reactions during lignin 
torrefaction, shifting the activation energy distribution of lignin to¬ 
ward lower energies [17], 

To formulate a model for predicting the extent of torrefaction of 
aspen and birch, it was assumed that their torrefaction is the result 
of independent torrefaction of their cellulose, hemicellulose and 
lignin content. Using the optimized value found for the DEAM of 
cellulose, xylan and lignin torrefaction, the contribution of each 
compound (c,- in Eq. (6)) was obtained so that the model best fits 
the TG derived V/V* profile versus time at 240 °C, 260 °C and 


280 °C. These values were used to predict the profile of V/V* versus 
time for the torrefied birch and aspen at 250 °C, 270 °C and 290 °C. 
The model is in agreement with the experimental data within the 
temperature range from 240 °C to 290 °C as seen in Figs. 8 and 9 
(r 2 = 0.97 for both aspen and birch torrefaction). The optimized val¬ 
ues of cellulose, xylan and lignin contribution (Cj) on the volatile re¬ 
lease of birch during its torrefaction are 0.584, 0.212 and 0.204, 
respectively. In the case of aspen wood, the contributions of the 
abovementioned compounds in the volatile release during torre¬ 
faction are 0.534, 0.235 and 0.231, respectively. Note that birch 
is composed of almost 58.9% of cellulose, 20.3% of hemicellulose 
and 20.8% of lignin [51] (values normalized on ash- and extrac- 
tives-free basis). The computed contributions of each compound 
(c,) on the volatile release of birch are close to the contributions 
of each compound in the chemical composition of birch wood. 
From this, it may be concluded that cellulose, xylan and lignin 
are decomposed independently during birch wood torrefaction 
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Fig. 10. TG derived profile of volatile release versus time for birch torrefaction and 
the results of its simulated profiles using the actual chemical composition of birch. 



Fig. 11. TG derived profile of volatile release versus time for aspen torrefaction and 
the results of its simulated profiles using the actual chemical composition of aspen. 

and the presence of alkali and alkaline metals mostly found in the 
lignin structure do not affect the rate of cellulose and hemicellu- 
lose torrefaction. As can be seen in Fig. 10, replacing q in Eq. (6) 
with the chemical composition of birch could fit well the experi¬ 
mentally obtained trends of V/V* versus time at temperatures from 
240-290 °C (r 2 = 0.97). However, in the case of aspen wood, the 
calculated values of q for cellulose, xylan and lignin were respec¬ 
tively, 4%, 15% and 8%, and different from their contributions in 
the chemical structure of aspen (51.2% of cellulose, 27.7% of hemi- 
cellulose and 21.1% of lignin [52]). Replacing q with the chemical 
composition of aspen leads to slight under-prediction of ver¬ 
sus time as shown in Fig. 11 (r 2 = 0.92). The difference between 
the chemical composition and the values obtained by the kinetic 
evaluation could be due to the effect of extractives decomposition 
on the torrefaction weight loss in the interrogated temperature 
range and/or to errors in the experimental techniques for obtaining 
the chemical composition of wood samples. 


5. Conclusion 

The torrefaction of two woody biomass samples (aspen and 
birch) and their main constituents (cellulose, xylan (representative 
of hemicellulose) and lignin) performed in a fixed-bed reactor 
showed that torrefaction improves the energy density of aspen 
and birch by reducing their oxygen and hydrogen content with re¬ 
spect to carbon in the solid residue. The equilibrium moisture con¬ 
tent of birch and aspen and of their main constituents (cellulose, 
xylan and lignin) was reduced by about 30-40% due to torrefaction 
and was explained according to FT1R spectra to an ameliorated 
hydrophobicity upon torrefaction of the treated samples. Further¬ 
more, the combined distributed activation energy model with a 
three-parameter log-normal distribution function could fit the 
experimental volatile release of xylan and lignin versus time while 
in the case of cellulose it exhibits some under-prediction. The dis¬ 
tributed activation energy models developed for cellulose, xylan 
and lignin torrefaction could also be combined in three-parallel 
reaction meta-models for the kinetic description of aspen and birch 
torrefaction. The model could fit the experimental data of birch 
and aspen torrefaction in the mild and severe torrefaction 
conditions. 
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